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When Channel Bonding is Beneficial for 
Opportunistic Spectrum Access Networks 

Shaunak Joshi, Przemyslaw Pawelczak, Danijela Cabric, and John Villasenor 

Abstract 

Transmission over multiple frequency bands combined into one logical channel speeds up data 
transfer for wireless networks. On the other hand, the allocation of multiple channels to a single user 
decreases the probability of finding a free logical channel for new connections, which may result in 
a network-wide throughput loss. While this relationship has been studied experimentally, especially 
in the WLAN configuration, little is known on how to analytically model such phenomena. With the 
advent of Opportunistic Spectrum Access (OSA) networks, it is even more important to understand the 
circumstances in which it is beneficial to bond channels occupied by primary users with dynamic duty 
cycle patterns. In this paper we propose an analytical framework which allows the investigation of the 
average channel throughput at the medium access control layer for OSA networks with channel bonding 
enabled. We show that channel bonding is generally beneficial, though the extent of the benefits depend 
on the features of the OSA network, including OSA network size and the total number of channels 
available for bonding. In addition, we show that performance benefits can be realized by adaptively 
changing the number of bonded channels depending on network conditions. Finally, we evaluate channel 
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bonding considering physical layer constraints, i.e. throughput reduction compared to the theoretical 
throughput of a single virtual channel due to a transmission power limit for any bonding size. 

Index Terms 

Opportunistic spectrum access, medium access control, channel bonding, performance analysis. 

I. Introduction 

In wireless communications networks, bonding multiple frequency bands into one virtual channel 
improves the transmission speed for individual network users at the medium access control (MAC) layer, 
as the channel throughput is theoretically linearly dependent on channel bandwidth. Practical evaluation of 
channel bonding concepts has attracted substantial amount of attention, especially from the wireless local 
area networks (WLAN) research community O-Q. Furthermore channel bonding is already present in 
some wireless networking standards, including IEEE 802.1 In, IEEE 802.1 lac and IEEE 802.22. 

Intuitively, in networks utilizing channel bonding, each data flow will be transmitted faster as more 
channels are used for an individual connection. However, more occupied channels for one user translate 
into less data transmission opportunities for other network nodes [6]. While this relationship is obvious, 
its analytical aspects have been relatively unexplored. Moreover, to the best of our knowledge, analyt- 
ical models exploring this tradeoff and potential benefits of channel bonding MAC protocols do not 
exist. Surprisingly, the channel bonding concept has not been addressed analytically in the context of 
Opportunistic Spectrum Access (OSA) network behavior, specifically considering MAC layer features 
(particularly random access mode) and PU and SU traffic characteristics. As OSA networks utilize 
multiple, potentially non-contiguous channels, with primary users (PUs) randomly reappearing on their 
respective channels, it is important to understand the performance of channel bonding MAC protocols 
for OSA networks. 

The objective of the paper is to understand the conditions under which it is beneficial to bond multiple 
frequency channels for OSA networks. To accomplish this, we present a mathematical framework based 
on a Markov chain analysis that enables the performance of channel bonding to be compared against 
classical OSA MAC protocols that do not use channel bonding. The framework enables the investigation 
of average MAC level network throughput as a function of the number of OSA network users, the number 
of frequency channels, channel bonding order, traffic level of secondary users (SUs), PU activity, and the 
spectrum sensing method, among many other parameters. The system model in our work is founded upon 
the model of Q and historically upon the model of [81]. This has provided for an ease of comparison of 
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our results with earlier models of (non-bonded) multichannel MAC protocols. However, unique properties 
of our model (especially the allowance of channel bonding, which has not been considered in Q, (H) 
resulted in a fundamentally new formulation of the Markov chain used in calculating the performance 
metrics. 

In the literature related to channel bonding performance, one of the first studies of channel bonding is 
found in ||2), where MAC protocols with channel bonding were advocated for throughput increase. It has 
been shown experimentally that channel bonding results in a large benefit for wireless networks. In this 
study only a small-scale network, i.e. with less than four nodes, was considered, thus it was unknown if 
the conclusions hold for large network sizes. In [3] the authors present measurement results of channel 
bonding functionality in IEEE 802.1 In. They show through experiments that when transmitted power is 
fixed for WLANs with and without channel bonding, the network with bonded channels obtains smaller 
throughput than a network without this option enabled. To alleviate this problem authors propose an 
algorithm that adaptively merges the channels depending on the perceived link quality, where users with 
the same link quality are associated with the same access point. The algorithm has been evaluated only 
through experiments. In another study Q a WLAN system with channel width adaptation was developed 
utilizing OFDM. A fixed number of channels were used, thus the effect of channel bonding adaptation on 
the network performance was not investigated. Summarizing, the investigation of channel bonding MAC 
protocols was strictly limited to the experimental platforms and it is unknown how channel bonding MAC 
protocols perform in scenarios not considered in the above mentioned studies. 

A separate group of papers on channel bonding performance in an OSA context relates to early 
analytical studies. The benefit of channel bonding (in a broader context of frequency agility) has been 
evaluated analytically in |9]. This study assumed that the OSA network assigns channels on a pre- 
detected pool of frequency bands, and MAC protocols are abstracted and not considered. In IfTOl a 
spectrum allocation framework for OSA-based WLAN networks utilizing dynamically changing channel 
bonds was presented. First, a static and centralized channel allocation framework was proposed. Then a 
distributed MAC protocol utilizing a spectrum allocation scheme called B-Smart was discussed, and an 
analytical model based on a Markov chain analysis was used to assess the control channel throughput of 
B -Smart. The model used approximating expressions related to the probabilities of frame collision, frames 
being idle and successful frame transmission from IfTTl . which were applicable only to non-modified IEEE 
802.11 systems with a distributed coordination function. The work in [10] has been extended to WLAN 
networks in (5|, where an analytical framework to evaluate efficiency of channel bonding for non-OS A 
networks with multiple access points has been proposed. However, as in the case of IfTOl , the framework 
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was applicable to static frequency allocation only, where the specifics of the protocol that network nodes 
use for communication were not a focus of the approach. The tool to model the problem was based 
on a linear programming formulation from which linear programming relaxation and a greedy heuristic 
were developed. A similar approach to analyze channel bonding was presented in Ifl2ll . where analysis 
of an OSA ad-hoc network with static channel allocation for users cooperating with access points was 
analyzed. Non-contiguous channel allocation was allowed with more than one channel assigned to a 
single SU. In summary, all existing analytical tools to evaluate channel bonding are applicable to static 
frequency allocation and do not consider MAC protocol features. In the context of our work we need to 
mention [13], where a channel bonding feature for IEEE 802.22 networks was investigated. However, the 
proposed model only used approximations to reach a tractable solution, was applicable to static resource 
allocation networks based on OFDMA, and did not consider random access features of the MAC protocol. 

Finally, we need to mention experimental OSA platforms with channel bonding. One of the first 
networks utilizing channel bonding for OSA network was the WhiteFi platform ['14'], i.e. a fully oper- 
ational WLAN network working in the TV white spaces, which utilized varying size channels for data 
communications, i.e. 5, 10 and 20 MHz. In WhiteFi channel bonding was possible only for adjacent 
channels. Performance of the platform was evaluated via simulations and experimentation with a limited 
set of prototype devices. A description of some of the algorithms and hardware used in WhiteFi platform 
was also discussed in [15]. Other OSA network prototypes utilizing channel bonding can be found in ifToll . 
The so called Jello framework fused multiple non-contiguous bands into one virtual channel, which was 
more flexible in this respect than WhiteFi. Indeed, the authors have shown a benefit of channel bonding, 
however measured only in terms of disruption rate and residual spectrum use. Implemented in a GNU 
Radio, Jello was evaluated in a testbed consisting of only a few nodes where the results were compared 
against simulation. In addition, off-line experimental studies of channel allocation for OSA, based on 
real-life measurement data of PU occupancy, have been presented in ifTTI . Again, while this interesting 
study clearly showed how networks can benefit from accessing PU occupied channels, and how channel 
bonding benefits OSA networks, there are no analytical results that were presented therein. 

The rest of the paper is organized as follows. The system model is introduced in Section [III while the 
analytical model is presented in Section [Till Numerical results are presented in Section ITVl and the paper 
is concluded in Section [V] 
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II. System Model 

We assume a single hop OSA ad-hoc network composed of N nodes communicating with each other 
using multichannel MAC with a dedicated control channel (DCC), as discussed first in Q, (H, [18l. 
Each node is assumed to transmit a saturated flow of framed data, where each new frame is generated with 
probability p > 0. Parameter p, i.e. the connection request probability, governs the collision resolution 
strategy. A higher value of p results in a higher number of collisions on the control channel, while 
a lower p results in fewer collisions leading to a lower channel utilization. Furthermore, note that the 
assumption of traffic saturation for each node follows the classical assumption of many other works 
that have considered performance analysis of MAC protocols. As |[T9l Sec. III-B] notes "(■■■) saturation 
throughput is a major performance measure to evaluate MAC protocols". We refer for example to |[20l 
Sec. II-C], ||2T1 Sec. IV], and [22, Sec. 5], where such an assumption has been made while analyzing 
OSA MAC protocols. Note that by "saturated" traffic we mean traffic which is always sent by every 
node, however, at non-regular intervals governed by the control channel access probability p. 

Nodes contend for the channel access by transmitting a Request to Send (RTS) frame via DCC, 
requesting a connection to other user of the network. The connection request is successful when only 
one SU requests a connectiory. The connection is established when the intended receiver responds with 
Clear to Send (CTS) frame on DCC when it is not involved in data exchange with other users. At the 
event of collision all connection requests are lost and users must contend for the DCC access after a 
random amount of time. This approach closely mimics the operation of the S-ALOHA protocol. 

PU and SU transmission is slotted (SUs follow PUs' slot boundaries). The assumption on the slotted 
transmission of the SU and its synchrony with the PU is common in theoretical MAC analysis for OSA 
networks, see for example [23, Sec. II], (21 Fig. 2], (25] Fig. 1], |[26j Sec. III-A] [27j Fig. 1], J28l Sec. 
II- A]. Moreover, the assumption of PU/SU slot synchrony is due to conformity with our previous work 
on multichannel MAC protocols for OSA networks, i.e. [7], [18], where the same assumption has been 
made. This allows for an easy comparison of the results obtained in this paper with the previous results we 

'The reason for a single hop setup was due to the conformity and ease of comparison with the earlier studies on multichannel 
MAC protocols listed above. This makes it possible to easily compare analytical results on different OSA MAC protocol flavors 
with those available in this paper. In addition, many earlier studies on channel bonding in OSA networks, including prototype 
evaluation of channel bonding, considered single-hop domain as well, see e.g. |14|, [15|, |17 |. 

2 Note that the application of our studies go well beyond single hop networks, as they are also applicable to centralized 
(access point-based) networks. The conversion of our analysis to centralized networks can be performed easily by adapting the 
construction of connection arrangement. 
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have obtained on OS A multichannel MAC protocol performance. Furthermore, certain papers, e.g. 1291 , 
ll30l do assume PU/SU synchrony in the practical network scenarios (in the case of [29]: secondary GSM 
network use, and in the context of [30]: secondary IEEE 802.16 network use). Then, we emphasize that 
the PU/SU slot synchrony assumption allows for obtaining upper bounds on the throughput in comparison 
to slot-asynchronous interface, as it has been remarked in f3~Tl . Please note however, that our model can 
indeed be generalized to a non-slotted PU/SU transmission. However, this would require further analysis 
on channel access policies, like those performed in fl32l . |[33l . which are beyond the scope of this paper. 

Each slot is T seconds long, with T s < T seconds of the spectrum sensing time performed at the 
beginning of the time slot by OSA network users and T — T s is equal to a length of one RTS/CTS frame 
exchange. Frame lengths are geometrically distributed such that the average frame length is l/q(k) time 
slots, where q(k) = {C(T — T s ) / d}k(3(k) > 1 is the probability of a time slot being occupied by the 
frame transmission.! C is the theoretical throughput of a single physical channel, k is the number of 
physical channels bonded to form a virtual channel, (3(k) is the throughput reduction factor for a k- 
bonded virtual channel, where j3(l) = 1 and Mk > 1 < p[k) < 1, described in detail in Section Hl-Cl 
and d is the size of the frame in bits. As we assume that the control channel throughput (in b/s) is the 
same as the throughput of data channel, i.e. C, we can easily calculate the length of the control packet 
as dp = {C(T — T s )}/p. In other words, in our model data and control packets, although of different 
lengths, are sent at the same rate. 

A. Multichannel Channel Bonding MAC Protocol Abstraction 

In the present paper, the main difference with respect to the original design of the multichannel MAC 
with DCC (see again Q, (H, lfl"8l ) lies in assuming that each new connection will utilize K out of 
M = Mt — 1 channels, where My is the total number of channel^) at time t, where K is the maximum 
channel bond order. In the case of less than K physical channels available for a new connection, all 
free channels are used for a frame transmission. We assume a set of accessible PU channels (which 
still can be randomly occupied by PUs) is known to the whole network and assigned unique identifiers. 
The sender, through the RTS frame, communicates to the receiver the information on the ordered set 
of available channels seen from its perspective, just like in [ 10 , Sec. 4.1] (in the context of continuous 



3 In other words q(k) is the probability that during frame transmission next time slot will be empty. This will end transmission 
and force a sender to contend over a control channel for a new virtual channel for new frame transmission. 

4 In a network with a DCC, from the pool of all channels Mt, one channel is used for control data exchange only. 
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frequency blocks). Next the receiver, through the CTS frame, replies with the set of ordered available 
channels from its perspective. After a successful RTS/CTS exchange both the sender and the receiver 
switch to the first K channels that are common to both of therd_ 

We assume that each communicating pair is able to bond non-adjacent frequency bands, with no guard 
bands considered, similar to (9|. Thus, our framework serves as an upper bound on the performance 
of any channel bonding protocol considering physical layer constraints on the channel bonding process. 
Note that an analysis of the optimal channel bond resource allocation, considering guard bands, were 
considered recently, e.g., in |[35l . In the remainder of this paper we refer to this channel bonding MAC 
protocol as flexible channel bonding. 

B. Primary User and Detection Process 

Each PU channel of bandwidth W at time slot t is randomly occupied by the PU with probability 
q p , where q p is geometrically distributed. While not all PU traffic is memoryless, refer e.g. to studies 
of IT361 that prove such a statement, many studies reveal that the complete opposite is true. Various PU 
systems can indeed be correctly described by a memoryless process. For example, in (371, one of the most 
complete long-term studies of channel occupancy by various wireless systems, the geometric distribution 
is quite common, constituting almost 60% of PU traffic distributions for systems of interest, including 
DECT, GSM (900/1800) and 2.4GHz ISnE. 

Due to a single hop domain analysis we can reasonably assume for tractability that each OSA node 
performs detection individually and the decision on the PU state is the same for each node in the OSA 
network. Furthermore, each OSA node is equipped with a wideband spectrum sensor to obtain information 
about all PU channels during the sensing phased Decision on the PU state is prone to errors due to 
false alarm, pt, or mis-detection, 1 — p^. Assuming Rayleigh fading and Gaussian noise the respective 
expressions are given in lfl"8l Eq. (1)] and [18, Eq. (3)]. For those expressions we will give respective 

5 In this paper we use RTS/CTS terminology to preserve the naming convention of |7|, |8|. However, note that the purpose 
of the RTS/CTS transmission in the context of our work is different than for multi-hop networks. For example, in IEEE 802. 1 1 
networks (and their derivatives) RTS/CTS mechanism is used to resolve the hidden/exposed terminal problem [34|. 

6 Note that the studies of non-memoryless distributions' impact on multichannel MAC protocols were considered in our earlier 
studies. We refer to (7] Sec. 5.2.6] for further discussion. 

theoretically, channel bonding would be possible with a narrowband/single antenna sensor, but only under certain limiting 
conditions. That is, slow PU activity with fast channel sensing time, i.e. orders of magnitude faster than the frame sending rate 
of the SUs. This would potentially allow to sense all PU channels in a round robin fashion by the SU. 
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values for the PU channel bandwidth W and SNR 7 of the PU signal at the SU energy detector with 
detection threshold 0,, resulting in required pf and pd, when presenting the numerical results. We assume 
that the OSA network loses throughput only due to false alarms, and it can successfully deconstruct a 
frame on the arrival of a PU during a mis-detection. On the other hand, we assume that the detection of a 
PU on any of the physical channels constituting one virtual channel causes termination of an entire frame 
transmission. Although this is a very conservative strategy it simplifies the analysis significantly, while 
also serving as a lower bound on the obtained MAC throughput considering frame disruption resolution 
strategies. 

C. Physical Layer Assumptions 

In the analysis we consider two cases with respect to virtual channel throughput. In the first case we 
assume that throughput of a virtual channel with k bonds is exactly kC, as in O, i.e. Vfc (3(k) = 1. In 
the second case we assume, due to the fixed power that users can emit per virtual channel, the throughput 
for a /c-bonded channel is (3(k)kC, where V/c > 1 (3(k) < 1 and /3(1) = 1. Note that the function j3(k) 
can be defined arbitrarily V/c > 1, depending on what physical layer constraint is considered, e.g., it can 
mimic the exponential decay of throughput for virtual channels of increasing orders, or it can mimic 
step-wise throughput reduction with increasing k. 

D. Comparison with other Multichannel MAC Protocols 

For comparison in this study we also consider other multichannel MAC protocols. The first protocol 
is a classic multichannel MAC protocol analyzed first in 03 j, where bonding is not allowed, i.e. channel 
bonding protocol with K = 1. The second protocol, denoted as i^-only bonding multichannel MAC, 
is a less flexible version of the main channel bonding MAC protocol. That is, when a newly admitted 
connection sees fewer than K free channels, then this connection is blocked and the user contends again 
for the channel access. This protocol mimics the operation of networks which are able to bond only 
a fixed number of channels per newly arriving frame. The third protocol is an adaptive version of the 
considered channel bonding multichannel MAC protocol, that changes the order of channel bonding 
dynamically depending on the network and/or traffic conditions. Details of the protocol will be presented 
in Section IIV-FI In the paper we will study the flexible channel bonding MAC protocol via analysis, 
while i<C-only channel bonding and adaptive channel bonding will be considered via simulations. 
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III. Analytical Model 



To calculate network-wide average MAC layer network throughput we propose to use a Markov 
chain analysis frameworlj^. The model presented here extends the work of |7] to the case of channel 
bonding. Therefore we follow the same naming convention and use the same definition of connection 
arrangement and termination probabilities (as the system model in this work follows those of (7 ] and the 
earlier model of (H). However, the structure of the Markov chain, formulation of transition probabilities, 
introduction of a new function: connection preemption probability, and calculation of performance metrics 
are fundamentally new and unique to our work. 

The roadmap of this section is as follows. First, in Section IIII-AI we introduce the definition of a 
Markov chain used to analyze our proposed MAC protocol. Next, in Section IIII-BI and Section IIII-CI we 
present the definition of the connection arrangement and connection termination probability, respectively, 
while in Section IIII-DI we introduce the definition of the connection preemption probability. These three 
functions govern the way connections are established and terminated within the OSA network of interest. 
Section ITlI-EI derives the complete Markov chain transition probabilities (using the three above mentioned 
functions) which are later used in Section IIII-FI to derive the performance metrics of interest. 

A. Definition of a Markov Chain 

Assuming without loss of generality that N > M, i.e. that the OSA network needs to compete for 
limited resources, let x t = (xi it , ■ ■ ■ ,xx,t) denote a state of the Markov chain at time t, where state, 
< Xi t t < [M/i\ denotes the number of active connections occupying i physical channels at time t. 
Then, X t = x^, . . . ,Xj is the set of all possible states for the considered system at time t, and y is 
the size of X t . Furthermore, let 




Pr(£i,m =«i, • • • , x K ,t+i = a K 



xi,t = b\ 



■ ■ -,XK,t = bic)- 



(1) 



Then the steady state probability vector 7r ; 



Xi,...,X K 



is obtained from transition probabilities as 



1 ' ' ■ Pi(x ltt =Xl,.. .,x K ,t = x K ), 



(2) 



while X = lim^oo X t 



Use of the Markov chain as an analytical tool to understand the operation of MAC protocols is a standard approach in 
communications literature, e.g. see recent work of |38|, |39| analyzing (non-channel bonding) MAC protocols for OSA. Therefore, 
we follow the same approach here. 
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B. Connection Arrangement Probability 

The probability that j new connections at time t+1 are successfully admitted through a control channel 
while m pairs of users actively exchange data on the virtual channel at time t is defined as |7 ( Eq. (15)] 

(N - 2m)p(l - p) N ~ 2rn - j , j = 1, 

where, without loss of generality, the common control channel is assumed to be absent from the PU 
activity. If one wants to assume that DCC is also randomly occupied by the PU, then (0) can be easily 
replaced with [7, Eq. (15) and (16)A 



c0') 



1 - s {1) 



C. Connection Termination Probability 

The probability that j connections at time t+1 finish transmission out of m connections using A;-bonded 
virtual channels at time t, is defined as |7, Eq. (14)] 

™)q(kY{l-q(k)) m -3, m>j>0, 



(k) T {j) 



0, 



otherwise. 



(4) 



D. Connection Preemption Probability 

The probability that existing SU connections are preempted by PUs is defined as 



>(/.*) 



Y 

E 

i=i 



K c i:j a. 



j=l m=l n=l 

xofd- 



EEE(,(«») + EE 



Cjf + l / 



m=l n=l 



1 



,Af- 



(5) 

where q c = q v pd + (1 — <Zp)p/ Eq. (13)] is the state of the PU observed by the SU networ]<0, z is the 
total number of PU occupancies on physical channels at time t, f = M — Y^=x w i tn gi ven i n ©> 
is the number of PU occupancies on physical channels not used by SU pairs, r t+ i = [ri^+i, . . . , rx,t+i], 



Note that the application of our model to ad ahoc networks allows for the analysis of infrastructure-based networks as 
well. For example, in the simplest case one needs to change the definition of Sm to (jV — m)p(l — p~j N ~ m ~^ for j — 1. 
This would mean that every new connection request involves only a single user (instead of two), e.g. a single user requesting 
connection to an access point. Obviously, other definitions of Sm can be used in our work that would reflect the operation of 
infrastructure-based networks better. For example, one can use |40 Eq. (3)], i.e. definition of a PRACH-like control channel 
used in 3GPP-based networks. 

l0 Note again that pf results from a certain energy detection process for a required pd\ see Section Hl-Bl for details. 
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where r&^+i is the total number of preemptions for connections using fc-bonded virtual channels at time 
t + 1, and Y is the total number of combinations of PU occupancy across all physical channels that 

account for r t+1 SU preemptions. Following Q and (O, to simplify the notation in (f5]) we have omitted 

(i i) 

time index t such that x = x t and r = r t +\. To define lm,n we introduce 

Ai,i • • • Aik Aik+i 

; •-. ; ; (6) 

Ay,i • • • Ay,k Ay,k+i 
where Vj : 1 < j < K Aij is a matrix, the rows of which represent the distribution of PUs that 
can occupy channels to ensure SU preemptions, and A ?i a'+i is a matrix, the rows of which contain 
the distribution of PUs on / unoccupied channels, while still maintaining rj preemptions. Specifically, 

Vj : 1 < j < K + 1, 



Ai, 



(id) 
1,1 



(id) 



(id) 

Cj ,1 



(id) 



(V) 



,(id) 



where i^n indicates the number of PU occupancies on one of cij currently active j-bonded connections 
in order to ensure a total of r-i SU preemptions across all such active connections. Note that oa'+i is 
the current number of free channels and is not related to the number of bonded connections, aj, defined 
V j : 1 < j < K. The term cy represents the total number of combinations of arranging PUs on aj 
currently active j-bonded connections to ensure r preemptions for the fth row of A. 



E. Transition Probabilities 

First, we introduce the indicator function: I(rj) = 1 if i] is true and I(rj) = otherwise, where rj 
represents the logical condition under test. We further introduce supporting notation that is used in the 
derivation of the transition probabilities: Ok denotes the vector of length K containing only zeros; Si 
denotes a vector of length K containing zeros except with a 1 value at index i; A = Y2i=i a i> which is 
the sum of all the current connections; and a = min (/, K). The transition probabilities are grouped into 
four broad cases shown below. 

1) If 3li : bi — ai = 1 and Vj^i: bj — aj = 0: This case deals with the event in which there is only 
one new connection being made while there are no terminating connections. Then 

r f +D = ]p) r (0) S ^pVf K \ (8) 
i=i 
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It is noted that the event is reflected directly in the right superscripts of the ^Ty Z ^ functions while there 

(v) 

is exactly one connection generation seen in the superscript of S x . As there are no terminations due to 
preemption by PUs, hence the use of all the all-zeros vector 0^-. 

2) If Vi : hi — cij = 0: In the second case we explore the event in which there is no change in the 
number of connections within a time slot. There are two subcases as listed below. 

First, if / > 0, deals with the condition in which the next state is not a state in which connections 
occupy all available channels. The state in which all channels are occupied is referred to as an edge state. 
Then 



,(*+!) 



i=l 



(9) 



where Mr^ is given as ([8]) replacing Snt with Sj?\ i.e. 



(i) r .(*+ 1 ) 



A' 



TT (i) 2<0)£WpCf,0 K ) 



(10) 



(2)-(*+l) 



i=l 
K-l 



rr (.i) T (o)(K) T (i) 



i=i 



v qW p(/,0jc) 



0, 

r A 



(3U*+1) 



n « T (0) 5 (D 

(/-q,0k+<5 q ) 



i=l 



x P. 



x+5 Q 



0. 



if / > K, 

z = 0, 

ax > 0, 
otherwise, 

z > 0, 
otherwise. 



(11) 



(12) 



Equation (fTOl corresponds to the event in which there are no terminating connections and no pre- 
emptions due to PUs. Equations (TTTT > and (fT2l are conditional and are dependent on whether PUs are 
absent or present, respectively, in order to preempt SUs. In (ITTb a SU bonding K channels terminates a 
connection while another SU generates a frame ensuring the same connection state. In (fT2l an additional 
SU with K-bonded virtual channel makes a connection, but is preempted by a PU. Note that the current 
connection state and the number of free channels is adjusted by a factor of a, which is the number of 
channels bonded where an additional connection is made. 
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Second, if / = 0, represents the transition necessary when the next state is an edge state. Then r t 
is defined as © with the respective probabilities defined as 

K 

(i) r (*+i) = Wrj 0) ( s^p^f^ +s ( x ) p^ f z ~ a,0l<) 



i=i 



(13) 



{2)Jt+l) 



K-l 



n(iW0)(Jf) T (l) 
-'-CM 1 a K 



8=1 



A 13 ^4 r X.2 J 



0, 



K 

if J] iftj < M, 

i=l 

M mod K = 0, 

ax > 0, 
otherwise, 



(14) 



and ( 3 )r^ t+1 ^ defined in the same way as (fl"2l ). Subcase Wr^* (defined as (fT3T )) represents the event 
in which there is no additional connections or terminations ensuring the same connection state, and also 
contains an excess condition where one additional connection is blocked because of all channels being 
occupied, i.e. the edge case. Subcase ( 2 V^ +1 ^ (defined as (fl4l) ) and subcase ( 3 )r^ <+1 ^ (defined in the same 
way as (fT2l) ) are conditional and represent the case in which an extra connection can be preempted by an 
SU and an additional connection with K bonded channels replaces a terminating connection, respectively. 

3) If 3i : bi — di < and z > 0: This case represents the event in which there is at least one 
terminating connection and preemptions are also possible because of the presence of PUs. 

Before we describe this case in detail, we introduce further notation. Let be a matrix where each 
row contains the combination of allowable terminations for the particular state transition, where the total 
number of such combinations equals c, i.e., 



© 



'1,1 



'1,2 



(15) 



"c,l "c,2 ■ ■ ■ "c,K 

where the rows of are denoted as [9\, . . . , 9k] T ■ Furthermore, let r = xt+i — Xj, where tj is the jth 
element of r. We consider the following three subcases. 

First, if / = and Vi, j : r 3 - — O^j > 0, there are terminations possible due to overlap, i.e. Vj : 
Tj — Oij > including terminations according to the rows of with both none and an additional excess 
connection generation, i.e. 



,(*+!) 



K 



i=l j=l 



(16) 
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Second, if / > and Vj : ta — 0{ j > and pi : bi 



a-, 



1, there are terminations possible due to 



3 v h3 

overlap and although free channels are available, an increase in connection does not occur. That is 

c \ K 

E II ) s { a p^-^r 1 j9i,j ' T ~ 

1=1 j=l 



,(t+l) 



i) 



K 



K 



K 



J=l 



j=i 



o(l)p(/+£f= 1 i*«.i-«.T-*0 Jf _, 

x (i) T (e,,i) . . .(«) . . . Wr(«*.«c) 



v(? (i) ,,(/•>:; - " 



6a) 



(17) 



In ( fTTl ) the logical argument in the outermost indicator function is used to determine if there are 
enough PUs present for preemption in the event that an additional connection generation occurs and 
terminations are insufficient. The logical argument in the innermost indicator function is used to determine 
if terminations for a-bonded connections alone are sufficient to reach the desired end state. 

And finally, if / > and 3H : 6j — a* = 1, there are free channels and there exists one and only one 
additional connection generation the following transition probability is used, i.e. 



„(*+!) 



£ 

i=i 



K 



K 



TT C?) r (*«,,) 5 (i)p(/+Ef =1 ^,„r-eo 

I 1 ttj ^4. x — Oi .z 



i=i 



(18) 



Note that as in (fTTT ) the indicator function is used to determine if preemptions combined with terminations 
are sufficient to reach the desired end state. 

4) If 3i : bi — < and z = 0: Similar to case described in Section HII-E3I this case covers the 
event in which a termination occurs, however there are no preemptions possible because there are no 
PUs occupying the system in the end state. That is 



r (*+i) _ (i) r (*+i) , (2) (m) 
r t — 't ^ r t 

For the purpose of simplifying the notation we introduce the following supporting function 



(19) 
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vP (/+Ef =1 ie i , i ,o K ) 



(20) 
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Then 



' t 



£(4 0) +4 1} 



if / > 0,aj = bj, 



5 



(i) 



^ 0, 



if / = 0, &j 
if hj bj , 
otherwise, 



(21a) 
(21b) 
(21c) 
(21d) 



where a, 



E. 7 =i 



j) and 6j = X^jLi The first two conditions, i.e. (12 lab and (12 lbb . represent 



the non-edge and edge cases, respectively, for when there are no additional generations in the system, i.e. 
the terminations subtracted from the current state is equal to the end state. The third condition, i.e. (12 lcb . 
represents the case in which there is an additional generation accompanying the one or more terminations 
in the current connection state. 
Furthermore 



(2U*+1) 



< 



(lM0«.i) . . . 



x {K) T (0i, K ) s m 



X P 



(M--£«jb 3 ,0 K ) 



if / > 0, 
$z:bi- t 



(22) 



1, 



x-0i,O ' 

0, otherwise. 
In (|22l there are additional free channels, but there is no increase in any of the bonded channels. In this 
case an additional connection generation in the a-bonded channel is terminated by adding a 5 a to the 
ath term of the termination vector row 6L 



F. Performance Metric Calculation 

Given the steady state matrix tt Xi ,...,x k we can compute the performance metrics of the system of 
interest. The system throughput is calculated as 



y K 



i=ij=i 

Note that the channel utilization can be calculated as 

r-r s E ie ^Ef=iXi,iivr(i) 



u 



(23) 



(24) 



T M 

Also note, that in d23l and (1241 both metrics are multiplied by the factor (T — T s )/T, due to incurred 
spectrum sensing overhead. 
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To conclude, in the context of our model the metrics of (1231) and (1241 ) are the fundamental descriptors 
of the channel bonding MAC performance in OSA. In the next section we will present numerical results 
to gain insight on the operation of such a protocol based on these two metrics. 

IV. Numerical Results 

We present a set of numerical results to evaluate channel bonding for OSA from a MAC layer 
perspective. First, in Section lTV-AI we assess the impact of PU activity level q p , introduced in Section ITl-B I 
on the system throughput. Secondly, in Section IIV-BI we investigate the effect of the ratio of channel 
pool size to the total number of OSA network users on the system throughput. In Section IIV-CI we 
investigate the effect of individual virtual channel throughput, kCf3(k), introduced in Section Hl-CI on 
total system throughput. Then, in Section IIV-DI we consider the impact of a virtual channel disruption 
strategy on system throughput, followed by studies on adaptive channel access control in Section IIV-EI 
In Section IIV-FI we discuss the construction of an adaptive channel bonding MAC protocol. Finally, 
in Section IIV-GI we present results on the effect of non-geometric PU duty cycle distribution on the 
performance of channel bonding. In the subsequent sections we focus on the throughput metric R solely, 
since channel utilization U is directly proportional to R. Since the issues of random access phase and the 
effect of spectrum sensing layer were investigated in earlier studies of multichannel MAC protocols Q , 
[18], we do not explore them here. 

For the sake of clarity and without loss of generality, we present the results by grouping them into 
two network scenarios unless otherwise stated: (a) PU channel pool M = 4, number of users N = 12, 
and frame size d = 5kB, representing a small-scale network; and (b) M = 12, N = 40, d = 20 kB, 
representing a large-scale network. The remaining parameters common to both scenarios are (unless 
otherwise stated): control channel access probability p = e~ 1 /N (as in Q, lfi~8l ), physical channel 
throughput C = 200 kB, slot length T = 1ms, sensing time T s = 0.1ms, throughput reduction function 
Vi f3(i) = 1 (perfect bonding scheme), 7 = OdB PU signal received power at each SU detector, detection 
threshold Q = 17.8 dB and W = 200 kHz bandwidth for which p d = 0.9 and p f = 0.02. We have tested 
the performance of channel bonding schemes for realistic bonding order, i.e. K G {2, 3}, as investigated 
in, e.g., @-a. 



A. Impact of PU Activity on System Throughput 

We evaluate the impact of PU activity q p on the system throughput R. The results are presented in 
Fig. Q] The channel bonding scheme considered in the analysis is represented in Fig. l(a)| and Fig. |l(b)| 
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(a) M = 4, N = 12, d = 5kB 



(b) M = 12, N = 40, d = 20 kB 



Fig. 1. Impact of PU activity level q p on system throughput for (a) small network and (b) large network; common parameters: 
p = e~ x /N, C = 200 kB, T = 1ms, T s = 0.1ms, Vz f3{i) = 1, 7 = OdB, n = 17.8 dB, W = 200kHz, p d = 0.9, p f = 0.02; 
N: i^-only bonding (simulations), A: flexible bonding (simulations), An: analysis. 



as 'A'. For comparison we present the results for the K-only channel bonding scheme (represented in 
Fig. |l(a)| and Fig. |l(b)| as 'N' and described in Section III-Db . where a new connection is dropped if and 
only if the number of available PU channels is less than K. 

We observe that the impact of PU activity level q p on system throughput R is strictly non-linear 
for every channel bonding order and network size. This behavior is due to the increased impact of 
PU preemptions on the active connections using virtual channels. Furthermore, increasing PU activity 
results in less throughput reduction. This is because for every channel bonding order just one active PU 
on any physical channel constituting a virtual channel causes disruption. Our result closely resemble 
simulation results presented in lfl4l Fig. 11], where with an increase in the number of access points 
generating background traffic on TV UHF channels the throughput of WhiteFi platform decreases semi- 
exponentially. Both Fig. 1(a) and Fig. |l(b)| present similar curve shapes. The largest impact is observed 
for low values of PU activity, i.e. < q p J$ 0.2. Most importantly, we see that for these two network 
scenarios, small- and large-scale, the benefit of channel bonding diminishes as the PU activity increases. 

For the small-scale network, the non-bonded system always outperforms the systems using channel 
bonding schemes. Furthermore, we observe that K-only bonding with K = 3 yields the lowest throughput, 
irrespective of the PU activity level, and this difference is largest for low q p values. Note that both the 
flexible and K-only schemes exhibit the same performance for K = 2 for the entire range of PU activity 
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because each scheme bonds the same number of channels regardless of how many of them are available. 
This is true for any scheme in which M is perfectly divisible by K. For other bond orders, flexible and 
K-only bonding schemes converge to the same value with an increase in q p . 

In contrast to the small-scale network scenario, we observe a benefit for channel bonding in comparison 
to a non-bonded system for < q p J$ 0.1 in the large-scale network scenario. But as q p increases we 
observe that the benefit of bonding diminishes just as in the small-scale scenario and eventually the 
system experiences a slight degradation in performance as compared to the non-bonded system. This 
degradation is attributed to an increase in PU activity, which leaves more channels unutilized due to 
more preemption of bonded connections by PUs. 

Finally, the analytical results obtained using the proposed mathematical model have been verified via 
simulations implemented in Matlab. We observe a perfect match between simulation and analytical results, 
as seen in Fig. |l(a)| and Fig. |l(b)[ where analysis is marked as 'An'. This validates the correctness of 
the analytical model. 



B. Impact of User Pool Size on System Throughput 

Exploring the performance of the protocol with a wider set of parameters, in contrary to Section ITV-AI 
we investigate the performance of the channel bonding protocol with two additional values for time slot 
length, i.e. T = 2ms (i.e. time slot length in HSDPA standard) and T = 5ms (i.e. time slot length 
in MobileWiMAX standard) |@T] Sec. IV-A], two various PU activity levels, q p G {0.05,0.1} (which 
are typical values for moderately used PU channels |[37l ), and two network sizes, small with M = 4 
and large with M = 8. Other parameters are the same as in Section IIV-AI We divide the results into 
two groups: Fig. |2(a)| and Fig. |2(b)| present results for a small-scale network and large PU activity, and 
Fig. 2(c) and Fig. |2(d)| present results for a larger-scale network and smaller PU activity. 

Interestingly, as the number of users in the system increase, certain channel bonding orders yield better 
results than others for a fixed PU activity level. This phenomenon is scenario-dependent. In the case of 
the small-scale network, when the number of users per channel is relatively small, i.e. N < 12, the system 
with K G {2,3} outperforms the system with no bonding, but only in the case of T = 5 ms time slot 
and small OSA user population size, compare Fig. 2(a) with Fig. |2(b)| As the number of users exceeds 
12, the throughput from the system without channel bonding begins to exceed that obtained from higher 
bond orders for both small and large time slot lengths, compare again Fig. |2(a) with Fig. |2(b)| This is 
attributed to the increased number of collisions during the random access phase that higher bond order 
systems experience as compared to the non-bonded system. This occurs because more SU pairs remain 
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(c) M = 8, T = 2 ms, q p = 0.05 



(d) M = 8, T = 5 ms, g p = 0.05 



Fig. 2. Impact of user pool size TV on system throughput for (a) (b) small-scale network and large PU activity and (b) twice 
the size of the small-scale network and smaller PU activity, using flexible channel boding scheme; set of common parameters 
is the same as in Fig. Q] assuming d = lkB; S: simulations, An: analysis. 



unconnected. In other words, a smaller number of connected SU pairs, each occupying more resources 
with higher bond orders, is worse off than more connected SU pairs, occupying fewer resources, for the 
non-bonded system. This phenomenon becomes more exaggerated as the number of users in the system 



increase. In the case of the large-scale network, see Fig. |2(c)| and Fig. |2(d)[ the non-bonded system results 
in a drastically reduced throughput in comparison to the other bonded systems. With a small number 
of users in the system and K = 3 bond order, throughput is maximized but again only for the case of 



T = 5 ms time slot, see Fig. 2(d) As the number of users increases, the system with lower bond orders 
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obtains higher throughput, as seen at N = 18 when K = 2 curve exceeds K = 3 (Fig. 2(c) with a 



small time slot length) and at N = 22 (Fig. |2(d)| with larger time slot length). The intersection of the 
throughput curves for the non-bonded system and the scheme with K = 2 occurs at a very large number 
of users because the rate in which throughput curves saturate is slower for networks with large channel 
pools. 

Our analytical tool demonstrates the importance of determining channel bonding performance, since 
the intersection between different channel bond orders is strictly dependent on the scenario and is very 
difficult to deduce intuitively. Again, as in Section IIV-A1 the analytical results have been verified via 
simulations. A very good match between simulations and analysis further confirms the accuracy of the 
proposed mathematical model. 

C. Impact of Virtual Channel Throughput on System Throughput 

We investigate the effect of virtual channel throughput on R as shown in Fig. [3l where we investigate R 
as a function of the SU frame size d. It has been suggested that the virtual channel throughput for higher 
bond orders is less than that of lower bond orders due to the increase in overhead when bonding [3]. We 
select three functions, as an example, to govern the way that virtual channel capacity is determined. These 
are: (a) perfect virtual channel throughput, for which VA: /3(k) = 1 (used in earlier two sections); (b) 
residual decrease in virtual channel throughput for which f3{h) = l/k a , a = 0.1; and (c) large decrease 
in virtual channel throughput for which (3(k) = l/k a , a = 0.5. The function /3(k) has been selected such 
that it decreases throughput exponentially and it penalizes channels with higher bond orders more than 
those with lower bond order^H. 

As the frame size becomes larger, the system throughput saturates. This observation is common for 
all considered scenarios, channel bond orders and virtual channel reduction factors. Note that the impact 
of frame length on system throughput in the context of non-opportunistic spectrum access multichannel 
MAC has been explained in much larger detail in [8, Sec. 4.4], with such discussion later extended 
to opportunistic spectrum access multichannel MAC in [7, Sec. 5.2.2] and [18, Sec. Hi-El]. We refer 
to these papers for a longer explanation of this phenomenon. To summarize these findings, with longer 
frames SUs occupy PU channels for a longer time on average resulting in an improved system throughput. 
Furthermore, as frames/packets get longer the SUs do not have to contend for resources as often via the 



"Note, again that any function /3(k) can be selected for evaluation. The particular function in this paper was considered to 
represent the typical wireless transmission scenario [3|. 
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Fig. 3. Impact of virtual channel throughput on system throughput for (a) small network and (b) large network for two channel 
bonding schemes K £ {2, 3} and three virtual channel throughput reduction factors, using flexible channel boding scheme: (a) 
perfect throughput (a = 0), (b) residual decrease (a = 0.1) and (c) large decrease (a = 0.5); set of common parameters is the 
same as in Fig. Q] assuming q p = 0.1. 
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Fig. 4. Impact of penalty factor on system throughput for (a) small network and (b) twice the size of the small network for 
two channel bonding schemes K £ {2, 3} and two values of PU activity q p £ {0.05, 0.1} as a function of penalty factor, using 
flexible channel boding scheme; set of common parameters is the same as in Fig. [2] with d = lkB and T = 2 ms as in Fig. |2(a)| 
and Fig. [20 



22 



control channel, therefore the number of collisions on the control channel become smaller, which again 
translates into improved system throughput. The increase in system throughput with an increase in frame 
size occurs despite the increase in collision probability between PU and SU transmissions (note that the 
per time slot probability of collision between PU and SU transmission is independent of SU frame/packet 
size). 

When frame sizes are relatively small, in the range of [10,20) kb, the impact of a residual penalty for 
virtual channels on system throughput does not drastically affect the performance as expected. However, 
with a larger penalty, system throughput performance greatly deteriorates, proving no additional benefit 
of channel bonding in an OS A network. The difference between system throughput with a = 0.1 and 
a = 0.5 is more profound for the small-scale than for the large-scale network scenario, compare Fig. |3(a) 



and Fig. |3(b)| Irrespective of the scenario, a bond order of K = 1 (i.e. non-bonded system) outperforms 
the other bonding schemes for all SU frame lengths for every penalty factor. Please refer to the earlier 
explanation in Section IIV-AI 

In another set of experiments, we investigate the effect of increasing the severity of the penalty a 
in f3(k) = l/k a on virtual channels by observing the overall throughput for different channel bond 
orders and two different PU activity levels, q p G {0.05, 0.1} with results presented in Fig. [4] We observe 
an approximately linear decrease in system throughput with an increase in the penalty factor a. For a 
small-scale network, see Fig. |4(a)[ as well as another network which is twice the size of the small-scale 
network, see Fig. |4(b)[ even a very small penalty for the virtual channel results in system throughput loss 
as compared to the non-bonded system, for example compare K = 1 and K = 2 for a ^ 0.02 and both 
values of q p . In the small-scale scenario, the throughput of higher bond orders is worse than that of the 
non-bonded system, as shown in Fig. |l(a)| However, in the case of the larger network, see Fig. 4(b) the 
non-linearities of virtual channel throughput are better accommodated when PU activity is low (again, in 
our case q p = 0.05). Hence the bonded systems provide for an improvement in throughput in the range 
of < a ;$ 0.04. This proves that there is a benefit of channel bonding even when the OSA network 
is unable to exploit the full theoretical capacity provided by the virtual channel. Finally, we add that 
(not shown in the Fig. |4] due to space constraints) as the PU activity increases, the effect of the penalty 
function becomes marginal as preemption dominates incoming SU connections. 



D. Impact of Virtual Channel Disruption Resolution Strategy on System Throughput 

Thus far we assume that once the PU is detected on any of the virtual channels the transmitted frame 
is lost. This approach serves as a lower bound on the channel bonding MAC protocol performance, refer 
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to Section III-BI for details. In this section we will relax this assumption and consider an example of a 
more flexible frame disruption resolution strategy. 

Specifically, we extend our simulation platform and consider a strategy denoted as Channel Switching 
(S), where the OSA network frame on the arrival of the PU (on any of the physical channels) is being 
switched to other free PU channelo- In other words, any frame transmission that uses i<C-bonded channels 
must be switched to another set of K unoccupied physical channels. Thus, the frame is lost only when 
not enough physical channels are found for switching. As S strategy is equivalent to B0S1 in Q we 
refer to (7J Theorem 1] proving that distributed channel switching can be performed by the OSA network 
without additional signaling transmission. Furthermore, (7J Sec. 4.3.2] describes in detail how the OSA 
multichannel MAC protocol with DCC resolves collisions when multiple sender/receiver pairs decide to 
switch to the same pool of physical channels. For notational convenience, in the remainder of this section 
the fundamental frame disruption resolution strategy (considered in the analysis) is denoted as A. 

In this section, in addition to throughput we investigate collision probability, /, defined as the probability 
of using a physical channel when the PU is present. Note that due to the channel switching process, 
to keep the comparison with flexible OSA channel bonding MAC fair, we virtually prolong the frame 
length by the channel switching delay T p such that q(k) = {C(T — T s + T p )/d}kf3(k). The results are 
presented in Fig. [5] where for fixed N = 2M, T p = 100 fis, and two values of q p and d we observe the 
performance metrics of interest. 

It is immediately observed that with the S strategy the OSA network with the channel switching 



performs much better than with the A strategy, see Fig. 5(a) and Fig. 5(c) (the increase is in orders of 
magnitude in comparison to strategy A). But on the other hand we also observe a large increase in the 
respective collision probabilities with the introduction of strategy S, see Fig. |5(b)| and Fig. |5(d)| This 



result clearly shows that there is a trade-off between frame disruption resolvability and collisions induced 
to PUs by OSA operation. The larger the OSA frame and and the larger the PU activity per channel, the 



larger the collision rate, see Fig. |5(d)| where almost 2% collision probability is reached (which is again 
in the orders of magnitude larger than in the A strategy). Interestingly, with increase in K the collision 
probability rate decreases. This is because with a large channel bonding pool, there is less probability 
that free channels can be found for switching and the frame is dropped, in-turn reducing probability of 

12 We emphasize that we consciously do not consider frame disruption resolution strategies which exclude physical channels 
occupied by PUs from the virtual channel, as this would result in an unfair comparison with the flexible case of the OSA channel 
bonding MAC protocol. This is because the frame transmission rate in the flexible case does not change dynamically, i.e. q(k) 
being a function of bonding size K is kept constant throughout the whole network operation. 
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Fig. 5. Impact of PU channel pool size on system throughput and collision probability for (a), (b) q p = 0.1, d = 5kB and (c), 
(d) q p — 0.3, d — 20 kB for two OS A frame disruption resolution strategies, fundamental (lower bound) (A) and switching (S), 
three different bonding levels K 6 {1, 2, 3}, N = 2M and channel switching time T p — 100 [is; set of common parameters is 
the same as in Fig.Q] 



further collisions with the PU. 

We want to emphasize that with the S strategy, the conclusions obtained for OSA channel bonding 
MAC with the A strategy still hold. Mainly, the order of bonding preference is preserved and in most 
cases K = 1 outperforms other bonding schemes in terms of throughput. 

E. Impact of Varying Channel Access Control and Prioritization on System Throughput and Fairness 
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1) Non-Uniform PU Channel Utilization: In previous sections we have assumed that the channels are 
used equiprobably by each PU. In this case, the OSA network channel selection strategy does not have 
an effect on the MAC protocol performance as each selected channel has the same probability of being 
disrupted by the PU. However, as the channels become non-uniformly used by PUs, the selection of a 
successive channel for each new bond might have a profound effect on the OSA network performance. 
Therefore in this section we investigate the effect of channel selection strategies on the performance of 
the channel bonding MAC protocol with a non-uniform distribution of PU activity over the PU channel 
set. 

Intelligent channel bonding strategies, assuming side-note information on PU channel properties were 
investigated in ifTTl Sec. 5.2] (denoted therein as frequency bundling), based on a large data set of 
measured PU channels. We will also consider such a system in the context of our work. Specifically, 
we will assume that the OSA network is aware of non-uniform PU channel usage and knows the first 
moment of PU channel distribution for every channel i E {1, • • • , M}. As an example, for mathematical 
tractability and without loss of generality, we assume that channel i is used by PU with probability 



where A q denotes the per channel PU activity imbalance factor. Note that A q = denotes equal channel 
utilization by the PU for each channel and A q > denotes an exponential increase in channel utilization 
with increasing i. Equation (|25T ) allows for a fair comparison of different channel selection strategies 
against uniform PU channel usage, as the mean of set {q p , • • • , qp M ^} for any A q is equal to the mean 
channel utilization in the case of uniform PU activity over all PU channels, i.e. q p . Note, however, that in 
(|25T ) q p must be selected such that it does not violate PU channel usage properties, i.e. Mi < q p < 1. 

In the simulation experiment we consider two channel selection strategies: random (R) channel selection 
and Least-used (L) channel selection. In case of the L channel selection strategy the OSA network, 
knowing how channels are used by the PU, always selects channels that are on average lowest used by 
the PU for each new bond (to minimize probability of collision with the PU and prolong transmission 
time). Apart from observing the OSA network throughput, for a considered distribution of channel usage 
by the PU, we also observe how uniformly channels were used by the SUs in the OSA network. For that 
we adopt a common fairness indicator denoted as Jain's index, defined as 1421 Sec. 3] 




2 ' 



(26) 
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We present the results in Fig. [6] The L channel selection strategy significantly outperforms the R 
channel selection in terms of obtained throughput, see Fig. |6(a)| and Fig. |6(c)[ The gain from using the 
L channel selection increases as A q becomes larger. This phenomenon holds irrespective of the channel 
bonding size K and is more evident as network size increases. An obvious result is that for the L 
channel selection strategy the fairness is much lower than that of the R strategy since the R strategy has 
almost a perfect fairness of 1, especially for the small network scenario (see Fig. |6(b)| and compare with 
Fig. |6(d)| ). However, a more interesting (non-obvious) result is that with an increase of bonding size K 
the fairness of the L strategy increases. This is due to more channels being used by the MAC protocol 
simulatneously with a larger bonding order. Furthermore, for K > 1 fairness decreases as A q increases. 
As the PU activity distribution becomes skewer, i.e. A q becomes large, the fairness of both the L and R 
strategies converge. This is especially visible in Fig. |6(d)| in the case of a large OSA network because 
most of the PU activity is concentrated over a small pool of channels, while in other cases the channels 
are sporadically used by PUs. Interestingly, varying the channel selection strategy does not change the 
conclusions obtained earlier, e.g. K = 1 reaches highest throughput followed by K = 2 and then K = 3. 

2) SU Channel Access Prioritization: In addition, we consider a fundamentally different channel 
bonding OSA MAC protocol, where priorities within the SU network are allowed. In particular, we 
consider a flexible channel bonding MAC, where each SU is able to send two types of frames: (i) regular 
and (ii) high priority. 

The priority-enabling flexible channel bonding MAC works as follows. Each newly generated SU frame 
is a high priority frame occurring with probability p^. An active, regular SU connection is halted and 
displaced on the arrival of a high priority SU connection when all channels are already occupied by 
active SU connections. A set of B € {0, ■ ■ ■ , M} halted connections in total is allowed, where B refers 
to the size of the buffer. Note that while in the real-life the most realistic setup for the OSA network is 
B = M (as buffering occurs by halting the existing connection of one SU transmitter/receiver pair by 
means of observing connection requests on the control channel [7, Sec. 4.3.3]) in numerical evaluation 
we will also consider the effect of B < M on the system throughput. Furthermore, for consistency with 
earlier results presented in this paper we do not consider buffering of SU connections on the event of 
preemption by PU connection, noting that the detailed consideration of multichannel (non-bonding) OSA 
MAC protocols with frame buffering, considering PU preemptions, has already been considered in (7J 
Sec. 4.3.3]. If more than one channel is freed from PU occupancy or SU transmission (of any kind), then 
one of the buffered connections resumes transmission until either the frame is successfully transmitted or 
preempted by the presence of a PU. The connections to be resumed are selected randomly from the pool 
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Fig. 6. Impact of PU channel usage imbalance on system throughput and fairness for (a), (b) small network and (c), (d) large 
network for two channel selection strategies, random (R) and least used (L), three different bonding levels K 6 {1,2,3} and 
q p — 0.2; set of common parameters is the same as in Fig.Q] 



in the buffer and per time slot one connection reconnect is possible (for consistency with the connection 
arrangement policy on the control channel). If during a single time slot new connection wants to access 
the channel pool and there is an existing connection to be resumed, the new connection gets priority 
over the buffered connection only if it is a high-priority connection. For regular connections, the buffered 
connection gets priority to connect to the channel. Furthermore, note that existing connection in the buffer 
does not contend again for channel access through the control channel. When one of the SU connections 
needs to be buffered, it is also selected randomly from the set of existing connections. In practical OSA 
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Fig. 7. Fairness of prioritization in flexible multichannel MAC as a function of prioritization rate Ph- (a) throughput, (b) low 
priority SU packet buffering rate; for three bonding orders, K £ {1,2,4} and two different buffer sizes, B £ {2,4}. Small 
network size is considered, q p = 0.05 while the rest of the common parameters are the same as listed for Fig.fJJ 



this can be done by, e.g., implementing a rule of buffering first the connections of a SU sender that have 
been most buffered so far (each SU can easily track which connections were buffered by observing the 
control channel, refer to J7J Theorem 1] by analogy). 

We consider fairness by observing simultaneously, for one (small scale) network seturo: (i) throughput 
obtained for the complete OSA network and (ii) average number of low priority SU frames (as a fraction 
of total SU frames) that are displaced to the buffer, both as a function of p^ The results are presented 
in Fig. |7] for three bonding orders, K G {1,2,4} and two different buffer sizes, B £ {2,4}. Observing 
the first metric, we see that the system throughput has a hyperbolic shape, maximized around p^ = 0.5 
and equal to system throughput with B = at p^ € {0, 1}, since it is at these two points that buffering 
is not effective. That is, when there are no high-priority connections (ph = 0) or when every connection 
is high-priority and cannot be displaced (p^ = 1). We observe that the buffer size B does not help in 
obtaining higher throughput for given value of K. Furthermore, an increase in K results in a decline in 
throughput due to a decrease in the bond adaptivity of each generated connection. 

The second metric, i.e. the fraction of total packets buffered presented in Fig. |7(b)[ is an important 
indicator of performance trade-off for the system. As the fraction of total buffered connections (B) 



l3 Note that we do not present accompanying results for large scale OSA network, as the buffering happens seldom in this 
case becoming statistically negligible. 
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increases, more connection disruptions occur for regular SU connections. Furthermore, we observe a trade- 
off with respect to ph since, for each curve, the highest fraction of buffered packets occurs approximately 
in the range of p^ G {0.4, 0.6}, causing SUs to occur maximum delay in reconnection. When observing 
frame buffering rate as a function of ph, we see that the considered MAC protocol is the most fair for 
Ph € {0, 1} demonstrating that the overall increase in system throughput incurs a penalty to regular 
connections. The shape of the packet buffering function is different for different combinations of B 
and K. That is, the smaller the K, the more skewed towards higher ph the function is. With higher K 
the effect of increasing B and its related increase in fraction of packet buffers becomes more profound 
(conversely, observe that curves for 'B = 4, K = V and 'B = 2, K = V completely overlap). 

The main message is to keep B < M for the same value of K, as the proposed MAC protocol obtains 
the same throughout for the same K irrespective of B (see again Fig. |7(a)[ ), while the packet buffering 
rate increases as B tends to M. Mean waiting times for the buffered connections are approximately equal 
for each case considered, specifically: (i) 4.59 slots for B = 4, K = 1, (ii) 4.55 slots for B = 2, K = 1, 
(ii) 4.73 slots for B = 4, K = 2, (iii) 4.49 slots for B = 2, K = 2, and (iv) 4.65 slots for B = 4, 
K = 4. 

F. Impact of Traffic/Network Dependent Adaptive Channel Bonding on System Throughput 

In the following study we extend the operation of the principal (flexible) channel bonding multichannel 
MAC protocol and K-on\y channel bonding multichannel MAC protocol with specialized channel adap- 
tation features. In other words, the bonding order K is adapted in an effort to maximize SU throughput. 
Such an observation can be made while observing Fig. |2j where depending on the user pool size, N, 
certain bonding orders obtain higher throughput than others. 

Let us assume that within a time interval i PU/SU traffic parameters are stationary. Then we can define 
the following optimization function for time interval i, necessary to find the the optimal bonding scheme: 

arg max R (27 a) 

K 

given M, N, d, p, q p , (27b) 

such that K < K mSLX ,p d > p d , r ,Pf <Pf,r, (27c) 

where K = Ki refer to the optimal bonding order at a specific value of q p for time interval i, K m3iX 
represents the maximum bonding level allowed, while p^ r and pd, r denote the OS A network- wide required 
probability of false alarm and mis-detection, respectively. Note that R is a function of all parameters 
considered in d27b| )- (|27c| ), while the framework to derive R is presented in Section JII] 
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To explore the effectiveness of the proposed optimization framework, we consider two network sce- 
narios identical to those considered in Section IIV-AI i.e., a small and a large network scenario, with 
-Kmax = 3, Pd, r = 0.9 and pf r = 0.1. All other parameters remain the same as in Section HV-A| except 
that PUs change their activity level in the course of the simulation from q p = in increments of A = 0.05 
to q p = 0.1. A small range of q p values were chosen in accordance with Fig. [J where it is observed that 
higher bond orders have better performance at q p is 0.1. 

Furthermore, as a bonding-dependent parameter we consider PU activity. During each of the above- 
mentioned time intervals the bonding order is assumed to remain constant. We assume that ( denotes the 
number of time intervals considered, while each K{ in the bonding scheme {«!,••■ , k^} is the bonding 
order for an associated value of q p . Note that we do not adapt K on SU traffic conditions as they are 
more difficult to estimate in real-life network operation than PU statistic^. For example, length of SU 
frames are not known a priori to the receiver (such information is not sent in the RTS frame in our 
system model), thus other nodes cannot estimate for how long a virtual channel will be occupied by a 
certain sender/receiver pair just by listening to a DCC. 



In Fig |8(a)| we present the SU throughput performance for ( = 3 and three different bonding schemes 
for a small network scenario and corresponding q p values: an optimal bonding scheme of {1, 1, 1}, 
obtained via solving d27a| )- d27bl , and for comparison {3, 2, 1} and {3, 3, 3}. We immediately notice that 
bonding, i.e. K > 1 for any value of q p , does not provide any gains and is in fact more harmful in terms 
of SU throughput. This is attributed to the bonding order being nearly equal to the channel capacity of the 
system, M, therefore many incoming SU connections that are granted access through the control channel 
are blocked. Furthermore, the collision resolution strategy used in the system makes it unfavorable for 
higher bond orders to maintain a connection because if a PU happens to occupy any single physical 
channel in a bonded virtual channel the entire frame is preempted, as described in Section III-BI 

In Fig |8(b)| we represent the SU throughput performance in a large network scenario, for the respective 
q p values, considering optimal bonding scheme of {3, 2, 1}, again obtained via solving (127al )- (127bl) . and 
for comparison {3, 3, 3} and {1, 1, 1} bonding schemes. In the large network scenario M > K, therefore 
a few bonded connections do not occupy the entire system bandwidth making it possible for multiple 
bonded connections to be present. This phenomenon allows the adaptive bonding scheme to utilize higher 
bond orders in time intervals in which q p ~ in order to maximize the usage of channels. This is clearly 
visible as the adaptive bonding schemes {3,2,1} and {3,3,3} outperform the single-bonded scheme 

14 For a recent discussion of efficient PU traffic estimation we refer e.g. to |43 | and the references therein. 
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SU packet length d (b) x1Q 4 SU packet length d (b) x1Q 4 

(a) M = 4, N = 12 (b) M = 12, N = 40 

Fig. 8. Impact of an adaptive channel bonding scheme on system throughput for (a) small network and (b) large network for 
three adaptive channel bonding schemes, i.e. K = {3, 2, 1}, K = {3, 3, 3} and K = {1, 1, 1} and two channel bonding MAC 
protocols, i.e. flexible (A) and A"-only (N), as a function of the SU frame length d for PU activity of q p £ {0,0.05,0.1}; set 
of common parameters is the same as in Fig.Q] 



{1,1,1}. Furthermore, {3,2,1} outperforms {3,3,3} as it is able to reduce the bonding order with an 
increase in PU activity level. 

In addition, as concluded in Section IIV-AI the A-only bonding scheme performs worse than the 
flexible channel bonding scheme for small-scale networks. This is because of the possibility of using a 
fixed number of channels, i.e. K in this particular case, for each connection is smaller then that of a 



large-scale network (notice the large difference in performance for the two protocol types in Fig. |8(a)| ). 
However, as the system size increases, the difference between the flexible and the A'-only channel bonding 
MAC diminishes because there is a higher probability A-order connections are made due to increased 
contentions on the control channel. In the case considered in this paper, the difference between the two 
protocol options is almost negligible, see Fig. |8(b)| 

G. Impact of Heavy-Tailed Distribution of PU Traffic on the Performance of Flexible Channel Bonding 
Protocol 

So far we have assumed that the PU traffic is distributed according to the geometric distribution. In 
this section we relax this assumption by considering the impact of a heavy-tailed distribution on the 
channel bonding MAC protocol performance, in the same way as it was considered in Q Sec. 5.2.6] for 
the non-bonding multichannel MAC protocols. In particular, we assess how the long-tail distribution of 
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Fig. 9. Impact of different PU "on" and "off" times distributions on OSA flexible channel bonding MAC protocol performance: 
(a) small network, and (b) large network, as given in Fig. Q] Legend symbols 'E', and 'L' denote geometric and log-normal 
distribution respectively, while the first and the second parameter of the two-letter symbol denote "off" and "on" time, respectively. 



the PU activity affects the performance of the flexible channel bonding MAC protocol. As an example 
of a long-tail distribution we use the log-normal distribution, as it was concluded in j37l Table 3 and 
4] that measured PU occupancies can be described by such a distribution in approximately 40% of the 
cases for GSM 900 uplink, GSM 1800 downlink, DECT and 2.4 GHz UNII channels. 

We present the protocol performance for four different combinations of "on" and "off" time distribu- 
tions, considering geometric and log-normal distribution, i.e. all possible combinations of "on" and "off" 
times obtained in 11371 Table 3 and 4]. They are denoted symbolically as: (i) EE, (ii) EL, (iii) LE and (iv) 
LL, where 'E' and 'L denote the geometric and the log-normal distribution, respectively, while the first 
and second position within the two-letter expression denote "off" and "on" times, respectively. Just like 
in E Sec. 5.2.6], the parameter of the log-normal distribution was selected such that its mean value was 
equal to qi = l/q p for the "on" time and qi = l — l/q p for the "off" time. Because the log-normal distribu- 
tion is continuous, it was rounded to the nearest integer, with the scale parameter a = y^log (vi/cf + l) 

and location parameter \i = log ^cf ^Jvi + cf V where q = 1/qi, v\ = (1 — qi)/qf is the mean and 
variance of the resulting discretized log-normal distribution. Note, as in (7J Sec. 5.2.6], that the variance 
of the resulting discretized log-normal distribution is equal to the variance of the geometric distribution 
for the same mean value. 

The results are presented for two types of networks denoted symbolically as small and large (see 
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Fig Q] for details), respectively in Fig 9(a) and Fig. 9(b) The values of PU activity, q p € {0.1,0.2,0.3}, 
represent low activity rates, which are of the most interest for the network designer. Furthermore, for the 
sake of clarity of the presentation we consider only two bonding orders, i.e. K € {2, 3}. First, we observe 
that the impact of log-normal PU channel occupancy distribution on the system throughput is not large 
as one might expect (for example for q p = 0.3 and the large network, Fig. |9(b)[ the difference between 
the EE and LE cases is less than 7%). If one wants to select the best combination of PU distribution 
that maximizes system throughput, it would be LE followed by LL. The impact of the non-geometric 



distribution on throughput is more visible for large networks than for smaller ones, compare Fig 9(a) 



with Fig. |9(b)| as in large networks there are more opportunities to exploit larger number of channels for 
transmission. All bonding orders are effected equally by the non-geometric process of the PU occupancy, 
thus all conclusions drawn for the channel bonding MAC protocol using the geometrically distributed 
PU traffic will also hold for other combinations of PU traffic. 



V. Conclusions 

In this paper we have developed a detailed analytical model for assessing system throughput of an 
ad-hoc Opportunistic Spectrum Access (OSA) network with channel bonding. Our analysis has led to 
a set of important conclusions. Firstly, we note that there is generally a benefit of channel bonding for 
OSA networks which is most prominent for low Primary User (PU) channel activity and/or large PU 
channel pools. On the other hand, in certain cases, channel bonding might result in significant throughput 
loss in comparison to a classical non-bonded system, i.e. when the PU activity is large or the users per 
channel ratio is very large. Secondly, to be able to exploit channel throughput fully, a medium access 
control protocol should adaptively change the channel bond order K depending on the network conditions 
(in particular the number of available resources and currently contending users). Thirdly, with certain 
physical layer constraints on maximum channel capacity, e.g. due to the limited transmission power per 
channel, irrespective of its bandwidth, the benefit of channel bonding for ad-hoc networks is easily lost. 
However, networks with a small users per channel ratio are still able to obtain higher system throughput 
than a regular non-bonded system, even when these constraints are present, provided, again, that the PU 
activity level is small. 

References 

[1] S. Joshi, P. Pawelczak, D. Cabric, and J. Villasenor, "Performance of channel bonding for opportunistic spectrum access 
networks," in Proc. IEEE GLOBECOM, Anaheim, CA, USA, Dec. 3-7, 2012. 



34 



[2] R. Chandra, R. Mahajan, T. Moscibroda, R. Raghavendra, and P. Bahl, "A case for adapting channel width in wireless 

networks," in Proc. ACM SIGCOMM, Seattle, WA, USA, Aug. 17-22, 2008. 
[3] M. Y. Arslan, K. Pelechrinis, I. Broustis, S. V. Krishnamurthy, S. Adepalli, and K. Papagiannaki, "Auto-configuration of 

IEEE 802.1 In WLANs," in Proc. ACM CoNEXT, Philadelphia, PA, USA, Nov. 30 - Dec. 3, 2010. 
[4] K. Tan, J. Fang, Y. Zhang, S. Chen, L. Shi, J. Zhang, and Y. Zhang, "Fine-grained channel access in wireless LAN," in 

Proc. ACM SIGCOMM, New Dehli, India, Aug. 30 - Sep. 3, 2010. 
[5] T. Moscibroda, R. Chandra, Y. Wu, S. Sengupta, P. Bahl, and Y. Yuan, "Load-aware spectrum distribution in wireless 

LANs," in Proc. IEEE ICNP, Orlando, FL, USA, Oct. 19-22, 2008. 
[6] "WLAN channel bonding: Causing greater problems than it solves," Texas Instruments, Tech. Rep. SPLY003, Sep. 2003. 
[7] J. Park, P. Pawelczak, and D. Cabric, "Performance of joint spectrum sensing and MAC algorithms for multichannel 

opportunistic spectrum access ad hoc networks," IEEE Trans. Mobile Comput., vol. 10, no. 7, pp. 1011-1027, Jul. 2011. 
[8] J. Mo, H.-S. W. So, and J. Walrand, "Comparison of multichannel MAC protocols," IEEE Trans. Mobile Comput., vol. 7, 

no. 1, pp. 50-65, Jan. 2008. 

[9] L. Cao, L. Yang, and H. Zheng, "The impact of frequency agility on dynamic spectrum sharing," in Proc. IEEE DySPAN, 
Singapore, Apr. 6-10, 2010. 

[10] Y. Yuan, P. Bahl, R. Chandra, T. Moscibroda, and Y. Wu, "Allocating dynamic time-spectrum blocks for cognitive radio 

networks," in Proc. ACM MobiHoc, Montreal, Quebec, Canada, Sep. 9-14, 2007. 
[11] G. Bianchi, "Performance analysis of of the IEEE 802.11 distributed coordination function," IEEE J. Sel. Areas Commun., 

vol. 18, no. 3, pp. 535-547, Mar. 2000. 
[12] T. Bansal, D. Li, and P. Sinh, "Fairness by sharing: Split channel allocation for cognitive radio networks," Ohio State 

University, Tech. Rep. Department of Computer Science and Engineering TR28, 2010. 
[13] J. Park, P. Pawelczak, P. Gr0nsund, and D. Cabric, "Analysis framework for opportunistic spectrum OFDMA and its 

application to the IEEE 802.22 standard," IEEE Trans. Veh. Technoi, vol. 61, no. 5, pp. 2271-2293, Jun. 2012. 
[14] P. Bahl, R. Chandra, T. Moscibroda, R. Murty, and M. Welsh, "White space networking with Wi-Fi like connectivity," in 

Proc. ACM SIGCOMM, Barcelona, Spain, EU, Aug. 17-21, 2009. 
[15] Y. Yuan, P. Bahl, R. Chandra, and P. A. Chou, "KNOWS: Cognitive radio networks over white spaces," in Proc. IEEE 

DySPAN, Dublin, Ireland, EU, Apr. 17-20, 2007. 
[16] L. Yang, W. Hou, L. Cao, B. Y. Zhao, and H. Zheng, "Supporting demanding wireless applications with frequency-agile 

radios," in Proc. NSDI, San Jose, CA, USA, Apr. 28-30, 2010. 
[17] V. Kone, L. Yang, X. Yang, B. Y. Zhao, and H. Zheng, "On the feasibility of effective opportunistic spectrum access," in 

Proc. ACM IMC, Melbourne, Australia, Nov. 1-3, 2010. 
[18] P. Pawelczak, S. Pollin, H.-S. W. So, A. Bahai, R. V. Prasad, and R. Hekmat, "Performance analysis of multichannel medium 

access control algorithms for opportunistic spectrum access," IEEE Trans. Veh. Technoi, vol. 58, no. 6, pp. 3014-3031, 

Jul. 2009. 

[19] S. C. Jha, U. Phuyal, M. M. Rashid, and V. K. Bhargava, "Design of OMC-MAC: An opportunistic multi-channel MAC 
with QoS provisioning for distributed cognitive radio networks," IEEE Trans. Wireless Commun., vol. 10, no. 10, pp. 
3414-3425, Oct. 2011. 

[20] N. Baldo, A. Asterjadhi, and M. Zorzi, "Dynamic spectrum access using a network coded cognitive control channel," IEEE 
Trans. Wireless Commun., vol. 9, no. 8, pp. 2575-2587, Aug. 2010. 



35 



[21] X. Zhang and H. Su, "CREAM-MAC: Cognitive radio-enAbled multi-channel mac protocol over dynamic spectrum access 
networks," IEEE J. Select. Topics Signal Processing, vol. 5, no. 1, pp. 110-123, Feb. 2011. 

[22] Y.-Y. Lin and K.-C. Chen, "Asynchronous dynamic spectrum access," IEEE Trans. Veh. Technol., vol. 61, no. 1, pp. 
222-236, Jan. 2012. 

[23] S. Wang, J. Zhang, and L. Tong, "A characterization of delay performance of cognitive medium access," IEEE Trans. 

Wireless Commun., vol. 11, no. 22, pp. 800-809, Feb. 2012. 
[24] Y.-C. Liang, Y. Zeng, E. C. Peh, and A. T. Hoang, "Sensing throughput tradeoff in cognitive radio networks," IEEE Trans. 

Wireless Commun., vol. 7, no. 4, pp. 1326-1337, Apr. 2008. 
[25] A. T. Hoang, Y.-C. Liang, D. T. C. Wong, Y. Zeng, and R. Zhang, "Opportunistic spectrum access for energy-constrained 

cognitive radios," IEEE Trans. Wireless Commun., vol. 8, no. 3, pp. 1206-1211, Mar. 2009. 
[26] L.-C. Wang, C.-W. Wang, and F. Adachi, "Load-balancing spectrum decision for cognitive radio networks," IEEE J. Sel. 

Areas Commun., vol. 29, no. 4, pp. 757-769, Apr. 2011. 
[27] H. T. Cheng and W. Zhuang, "Simple channel sensing order in cognitive radio networks," IEEE J. Sel. Areas Commun., 

vol. 29, no. 4, pp. 676-688, Apr. 2011. 
[28] K. Bian, J.-M. Park, and R. Chen, "Control channel establishment in cognitive radio networks using channel hopping," 

IEEE J. Sel. Areas Commun., vol. 29, no. 4, pp. 689-703, Apr. 2011. 
[29] P. Papadimitratos, S. Sankaranarayanan, and A. Mishra, "A bandwidth sharing approach to improve licensed spectrum 

utilization," IEEE Commun. Mag., vol. 43, no. 12, pp. S10-S14, Dec. 2005. 
[30] P. Gr0nsund, H. N. Pham, and P. E. Engelstad, "Towards dynamic spectrum access in primary OFDMA systems," in Proc. 

IEEE PIMRC, Tokyo, Japan, Sep. 13-16, 2009. 
[31] J. Gambini, O. Simeone, Y. Bar-Ness, U. Spagnolini, and T. Yu, "Packet-wise vertical handover for unlicensed multi- 
standard spectrum access with cognitive radios," IEEE Trans. Wireless Commun., vol. 7, no. 12, pp. 5172-5176, Dec. 

2008. 

[32] S. Geirhofer, L. Tong, and B. M. Sandler, "Dynamic spectrum access in the time domain: Modeling and exploiting white 

space," IEEE Commun. Mag., vol. 45, no. 5, pp. 66-72, May 2007. 
[33] S. Huang, X. Liu, and Z. Ding, "Optimal sensing-transmission structure for dynamic spectrum access," in Proc. IEEE 

INFOCOM, Rio De Janeiro, Brazil, Apr. 19-25, 2009. 
[34] Z. J. Haas and J. Deng, "Dual busy tone multiple access (DBTMA) — a multiple access control scheme for ad hoc network," 

IEEE Trans. Commun., vol. 50, no. 6, pp. 975-985, Jun. 2002. 
[35] H. B. Salameh, M. Krunz, and D. Manzi, "Design and evaluation of an efficient guard-band-aware multi-channel spectrum 

sharing mechanism," University of Arizona, Tucson, AZ, USA, Tech. Rep. TR-UA-ECE-201 1-1, 2011. 
[36] S. Geirhofer, L. Tong, and B. M. Sadler, "Dynamic spectrum access in the time domain: Modeling and exploiting white 

spaces," IEEE Commun. Mag., vol. 45, no. 5, pp. 66-72, May 2007. 
[37] M. Wellens, J. Riihijarvi, and P. Mahonen, "Empirical time and frequency domain models of spectrum use," Elsevier 

Physical Communication Journal, vol. 2, no. 1-2, pp. 10-32, Mar. -Jun. 2009. 
[38] X. Y. Wang, A. Wong, and P.-H. Ho, "Stochastic medium access for cognitive radio ad hoc network," IEEE J. Sel. Areas 

Commun., vol. 29, no. 4, pp. 770-783, Apr. 2011. 
[39] P. Wang, D. Niyato, and H. Jiang, "Voice-service capacity analysis for cognitive radio networks," IEEE Trans. Veh. Technol, 

vol. 59, no. 4, pp. 1779-1790, May 2010. 



36 



[40] P. Zhou, H. Hu, H. Wang, and H.-H. Chen, "An efficient random access scheme for OFDMA systems with implicit message 

transmission," IEEE Trans. Wireless Commun., vol. 7, no. 7, pp. 2790-2797, Jul. 2008. 
[41] V. H. G. E. 0ien and D. Gesber, "Throughput guarantees for wireless networks with opportunistic scheduling: A comparative 

study," IEEE Trans. Wireless Commun., vol. 6, no. 12, pp. 4215-4220, Dec. 2007. 
[42] R. K. Jain, D.-M. W. Chiu, and W. R. Have, "A quantitative measure of fairness and discrimination for resource allocation 

in shared computer systems," Digital Equipment Corporation, Hudson, MA, Tech. Rep. DEC-TR-301, Sep. 1984. 
[43] W. Gabran, C.-H. Liu, P. Pawelczak, and D. Cabric. (2012) Primary user traffic estimation for dynamic spectrum access. 

[Online]. Available: [http://arxiv.org/abs/1203.2393l 



